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Abstract 



When perfectly conducting or insulating inclusions are closely located, stress which 
is the gradient of the solution to the conductivity equation can be arbitrarily large as 
the distance between two inclusions tends to zero. It is important to precisely char- 
acterize the blow-up of the gradient. In this paper we show that the blow-up of the 
gradient can be characterized by a singular function defined by the single layer potcn- 
Oh ■ tial of an eigenfunction corresponding to the eigenvalue 1/2 of a Neumann-Poincare 

type operator defined on the boundaries of the inclusions. By comparing the singu- 
■ lar function with the one corresponding to two disks osculating to the inclusions, we 

quantitatively characterize the blow-up of the gradient in terms of explicit functions. 
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^ '■ 1 Introduction 

(N 

Let Di and D 2 be bounded simply connected domains in Mr, d = 2, 3, whose boundary 
O ! regularity will be specified later. Suppose that they are conductors, whose conductivity 

CN ' isfc,0</c7^1<oo, embedded in the background with conductivity 1. Let a denote the 

conductivity distribution, i.e., 

a = kx(D 1 UD 2 )+ X (^ d \(D 1 UD 2 )), (1.1) 
?—i ■ 
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where x is t ne characteristic function. We consider the following elliptic problem: for a 
given entire harmonic function h in R rf , 



V-crVu = in IT, , . 

S rf (1-2) 

|u(x) - /i(x) = Odxl 1 -^) as|x|^oo. 

Let 

e := dist(Di, D 3 ), (1.3) 

and assume that e is small. We emphasize that the shapes of D\ and D2 do not depend 
on e. More precisely, there are fixed domains D\ and D2 such that Dj is a translate of 
Dj, namely, there are vectors ai and Si2 such that 

Dj = Dj + a.j, 3 = 1,2. (1.4) 

The problem is to estimate |Vit| in terms of e when e tends to 0, or to characterize the 
asymptotic singular behavior of Vu as e — > 0. 

If k stays away from and 00, i.e., c\ < k < C2 for some positive constants c\ and C2, 
then |Vii| is bounded regardless of e as was proved in [101 [T9l I18j . In fact, it is proved 
that the C 1,a norm of u is bounded regardless of e when dD\ and dD2 are C 2,ot smooth. 
However, if k degenerates to either 00 (perfectly conducting case) or (insulating case), 
the ellipticity holds only outside D\ and D2 and completely different phenomena occur. 

When k = 00, the problem becomes 



An = in R a \ (D x L>D 2 ), 

u = Aj (constant) on dDi, i= 1,2, (1-5) 
u(x) -h(x) = OGxl 1 -^ as |x| 



00, 



f du 




f du 




+ 


JdD 2 dvW 



where the constants Aj are determined by the conditions 

= 0, (1.6) 

with z/( J ) being the outward unit normal to dDj, j = 1, 2. In that case, Vu may blow up 
as e tends to 0. 

As shown in [T71 [TTJ El [23 [Ml E] , in two dimensions the generic rate of gradient blow- 
up is e -1 / 2 , while it is |elog e| _1 in three dimensions [3[8l[20]. The blow-up of the gradient 
may or may not occur depending on the background potential (the harmonic function h 
in (|1.2p ) and those background potentials which actually make the gradient blow up are 
characterized in [4] when D\ and D2 are disks. In two dimensions, the perfectly insulating 
case, where k = 0, can be dealt with using the conjugate relation (see [T71 [5]) and in this 
case the blow-up rate is also e -1 / 2 . It is a challenging open problem to clarify whether 
|Vn| may blow up or not in the insulating case in three dimensions and to find the blow-up 
rate if the blow-up occurs. It is also a quite interesting problem to clarify the dependence 
of |Vii| on k as k —> 00 or k — >■ 0. In this relation we mention that a precise dependence 
on k when D\ and D2 are disks was shown in [5l [3] . It is worth mentioning that a similar 
blow-up phenomenon for the p-Laplacian equation was investigated in 
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Above mentioned results for k = oo are estimates of |Vu| in upper and lower bounds 
sense, namely, 

Ci^(e) < ||Vu||oo < C 2 V(e) (1.7) 
for some implicit constants C\ and C 2 where ip(e) = e" 1 / 2 in two dimensions and ^(e) = 
|e log e| — 1 in three dimensions. The constant C\ can be zero or positive depending on the 
background potential h. In order to have a better understanding of the concentration of 
the gradient it is desirable to pursue deeper investigation on the blow-up nature of Vtt. In 
this direction there is a recent work [13] where the blow-up nature of Vu is characterized 
in terms of an explicit singular function. It is shown that if D\ = B\ and D2 = B2 are 
circular inclusions of radius r\ and r 2 , respectively, and k = 00, then 

u(x) = — — (n- V/i)(p)(ln x-pi - In x-p 2 )+r(x), (1.8) 
r\ + r 2 

for x € R 2 \ (B\ UB2), where pi £ D\ is the fixed point of R1R2 where Rj is the reflection 
with respect to dBj, j = 1,2, P2 £ B2 is the fixed point of R2R1, n is the unit vector in 
the direction of p 2 — Pi, and p is the middle point of the shortest line segment connecting 
dBi and dB 2 . In (|1.8p . Vr is bounded independently of e and thus the blow-up of Vn is 
completely characterized by the singular function 

?b(x) := — (ln|x-pi| -ln|x-p 2 |). (1.9) 

The purpose of this paper is to establish a characterization of the blow-up of Vu similar 
to (jl.8|) when D\ and D 2 are strictly convex simply connected domains in K 2 . In doing so, 
the Neumann- Poincare (NP) operator denoted by K* and defined on L 2 {dD\) x L 2 (<9D 2 ) 
plays a crucial role. The NP operator is a classical notion and appears naturally when 
we solve the boundary value problems using layer potentials. It also appears naturally 
when we solve the transmission problem (jl.2p . See the next section for a definition of the 
NP operator for the problem of this paper. This operator has 1/2 as an eigenvalue with 
multiplicity 2. (If the inclusion has N simply connected components, then the multiplicity 
of 1/2 is N .) If two inclusions are disks, then (J^fylaDi, ^(f) \dD 2 ) T ■> where qs is given by 
(|1.9p . is an eigenfunction of the NP operator on L 2 (dD\) x L 2 (<9D 2 ) corresponding to 1/2. 
(This fact was also observed in [9].) Here, T denotes the transpose. 

Let g = (g^\ g^) T be the eigenfunction of K* on L 2 (dD\) x L 2 (<9D 2 ) corresponding 
to the eigenvalue 1/2 and satisfying 

# (1 W = 1, / c/ 2 W = -l. (1.10) 
dD 1 JdD 2 

We will prove that such an eigenfunction does exist. Let h be the background harmonic 
function introduced in (11.21) and let 



h:= (1.11) 

h\dD 2 \ 

Let (h, g) be the inner product on L 2 (dDi) x L 2 (dD 2 ), i.e., 

(h,g)=/ hg W da+ [ hg^da. (1.12) 
JdD 1 JdD 2 

With these notions in hand, we can state the main result of this paper. 



3 



Theorem 1.1 Let D\ and D 2 be strictly convex simply connected domains in M. 2 with 
Q2,a smoo if l boundaries for some a G (0,1]. Let zi G clDi and Z2 G <9-D 2 be the closest 
points, and let e := dist(D\, D 2 ) = \z\ — z 2 |, K j ^ e ^ e curvature of dDj at zj, Bj be the 
disk osculating to Dj at zj, j = 1,2, and qs be the singular function in tl.9\) associated 
with disks B\ and B2. Then, the solution u to il.5\) satisfies 

«(x) = ~^lM=a e (Zfl(x) + r(x) J x G R 2 \ (D 1 U D 2 ), (1.13) 
■y/e(/ci + k 2 J 

where a e is a constant bounded independently of e of the form 

«. = { ItSfA' n ^ £( ,°' 1) ' — »• ( U4 > 

[ 1 + 0(|Veme|) */ a = 1, 

and 

l|Vr|| LOO{f7) < C (1.15) 

/or some constant C independent of e. Here Q = Q\ \ (D\ U D 2 ) and Vt\ is an open set 
containing D\ U D2 ■ 

We emphasize that (|1.13j) is a pointwise relation and hence describes the behavior 
of Vu(x) in terms of the gradient of the function g^(x). One can see from the explicit 
expression (|1.9p that |Vgs| attains its maximum at z\ and z 2 , and that 



|Vg B (z,)| = ^^^ + 0(l). (1-16) 

(See (|5.5p .) So, (j 1 . 1 3 1) shows that |Vtt| is bounded regardless of e if (h, g) = 0. Moreover, 
it yields a new improved estimate: 

llVulUoCR^fliui*.)) = ° e| Y +0(1), ase^O. (1.17) 

Since 

I (h,g) I <CVi (1.18) 

for some constant C independent of e as shown in (|5.24p . we can also infer from (|1.17D 
that the generic rate of blow-up is e -1 / 2 . 

The (global) strict convexity assumption of D\ and D 2 in Theorem 11.11 can be relaxed 
a little. Instead, if we let z = Zl jj" Z2 , then it is enough to assume that there are 5 and Co 
such that dDj n B$(z) is strictly convex for j = 1,2, and 

dist^Di \ B s (z), 8D 2 \ B s (z)) > cq. (1.19) 

This can be shown by exactly the same proofs as in this paper. 

We also obtain similar results for the insulating case and boundary value problems. 
The problem for the insulating case, obtained by taking the limit as k — > of (|1.2p . is 
given by 

An = in R d \ (D l UD 2 ), 

= on 8Di, i = 1,2, (1.20) 

+ 



n(x) - h(x) = 0(|x| 1_d ) as |x| -)■ 00. 



4 



If u is the solution to (|1.20p . then its harmonic conjugate u 1 - is the solution to (|1.5p with 
h replaced with its harmonic conjugate h^. Therefore, we may apply Theorem 11.11 to u 1 - 
to obtain an analogous result for u. 

On the other hand, if we consider the boundary value problem 

V-aVu = in ft, (1.21) 

with either Dirichlet or Neumann boundary conditions on the boundary 9ft of a smooth 
domain ft containing D\ and D2 and k — > or 00, then, in this case, the harmonic function 
h defined by 

h(x) = ~[ ln|x-y|^(y)^(y)-i- / (x ~ y ' ^ y)) u(y) da(y), (1.22) 

with v being the outward normal to 9ft, plays the role of h in the whole space problem 
(jl.5p or (|1.20p . A similar result on the characterization of the gradient blow-up can be 
obtained using exactly the same arguments as for the whole space problem. 

The results of this paper can be applied for solving two longstanding problems. The 
first one is the study of material failure. In fact, the problem of estimation of the gradient 
blow-up was raised by Babuska in relation to the study of material failure of composites 
[6]. In composites which consist of inclusions and the matrix, some inclusions may be 
closely located and stress occurs in between them. The problems fjl .2[) . (jl.5p and (|1.20p are 
conductivity or anti-plane elasticity equation, and Vu represents the shear stress tensor. So 
results of this paper provide clear quantitative understanding of the stress concentration, 
which will be a fundamental ingredient in the study of material failure. 

The second application is computation of the electrical field in the presence of closely 
located inclusions with extreme conductivities (0 or 00) which is known to be a hard 
problem. Because |Vit| becomes arbitrarily large, we need fine meshes to compute Vu 
numerically. Since (| 1 . 1 3 1) for example provides complete characterization of the singular 
behavior of Vu, the complexity of computation can be greatly reduced by removing the 
singular term there. In fact, effectiveness of this scheme is already demonstrated in |13j 
when inclusions are disks using (|1.8p . It is worth mentioning that unlike (|1.8p where the 
coefficient of qs is explicitly determined by h, computation of the constant a e (h, g) in 
(|1.13p may cause a problem when inclusions are of general shape. We will investigate this 
issue in a forthcoming work. 

This paper is organized as follows. In the next section we introduce the single layer 
potential and define the Neumann-Poincare operator. In section 3, we construct eigenfunc- 
tions of the NP operator corresponding to the eigenvalue 1 /2 and prove that its multiplicity 
is 2. In section 4, we construct a singular function using eigenfunctions constructed in the 
previous section and characterize the gradient blow-up in terms of the singular function. 
In section 5, we estimate the potential difference of the solution to (|1.5p . Section 6 is to 
prove Theorem 11.11 Sections 7 and 8 are for the insulating case and the boundary value 
problem, respectively. In the last section we prove a lemma used in Section 6. 
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2 Preliminaries 

Let D be a bounded simply connected domain in M. d , d = 2, 3, with a Lipschitz boundary. 
The single layer potential 5dM of a function ip £ L 2 (dD) is defined as 

5 D M(x) = / r(x - yMy) da(y), x e 

J9D 



where T is the fundamental solution to the Laplacian, i.e. 



r(x) 



Then, for ip E L 2 (dD), we have 



— In Ixl 

2vr 1 1 



47T 



d = 2 
cZ = 3 



(2.1) 



d_ 

du 



Sd[<p] 



x 



where 



sr(x 



dD 



dv(x) 



-^(y)da(y), x£^. 



(2.2) 



(2.3) 



Here, denotes the normal derivative and the subscripts + and — represent the limits 
from outside and inside D, respectively. The operator IC* D is called the Neumann-Poincare 
(NP) operator. 

We now consider the configuration where there are two inclusions D\ and D 2 which 
are closely located. Suppose that the conductivity of the inclusions is k 7^ 1 while that of 
the background is 1, so that the conductivity distribution is given by (jl.ip . For a given 
entire harmonic function h in R , we consider the problem (|1.2p . 

It is known (see for example \H\ [T5] ) that the solution u to (jl.5p can be represented 

as 

u(x) =/i(x)+5 Bl [^](x)+5,) 2 [^](x), x G R d \ (Di U D 2 ) (2.4) 

for a pair of functions (<p^\ <p^) £ Lq(ODi) x Lq(8D2) (Lq denotes the set of I? functions 
with mean zero). Since u is constant on dD\ and <9Z?2> we have 







dud) 



(S Dl [^]+S D2 [p i2) }) 



dh 
"dvW 



on dDj, j = 1,2, 



which, according to (|2.2p . may be written as 



5/1 



on dDi, 
on <9Z?2- 



(2.5) 
(2.6) 



Here, -^jj denotes the outward normal derivative on dDj, j = 1, 2. This system of integral 
equations can be written in a condensed form as 



1 



(-I-K*)[<p] = dh, 



(2.7) 
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where 



(with / being the identity operator), and 



I 


0" 







/ 


, IC* : = 



JC 



^D 2 



(2.8) 



dh 



dh 
dh 



If there are N simply connected inclusions, say D±, . 
NP operator K* is defined by 



Z?at, then the corresponding 



IC* 

a 



IC 



T)Sd 2 



8vW 



Sdi 



a 



dv( N ) 



D 2 



Sd 2 



L>n 



(2.9) 



We make note of some important properties of the NP operator IC* whose proofs can 
be found in [I]. Let U := L 2 (dL>i) x L 2 (dD 2 ) and n := Lq(<9£>i) x L§(9D 2 ). We have 



• IC* maps H into itself, and T~Lq into itself. 

• For any A with A < -1/2 or A > 1/2, AI - 

• ^1 — JC* is invertible on "Ho- 



is invertible on %. 



• All the eigenvalues of IC* belong to (-1/2, 1/2]. 

One fact of crucial importance is that IC* can be symmetrized. To see this we introduce 
the operator § acting on T~L as 

g = S Dl S D2 

It is worth making a remark on the operators off diagonal. For example, Sd 2 on the top 
right corner is an operator from L 2 (dD2) into L 2 (dD\). It is proved in pQ based on a 
result in [16] that — § is positive semi-definite and self-adjoint, SIC* is self-adjoint, and 
hence there is a self-adjoint operator A on % such that 



(2.10) 



V^SIC* = AV^S. 
In other words, IC* is self-adjoint with respect to the inner product 

(ip, V>)a := -(§[(p],tp), ip,ip en. 



(2.111 



(2.12) 



3 Eigenfunctions of K* 

We now construct eigenfunctions of IC* corresponding to 1/2. Our construction plays an 
essential role in understanding the blow-up of the gradient. We first prove the following 
lemma. 
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Lemma 3.1 For i = 1,2, there is a unique solution vi to 



' Avi = inR d \D 1 l)D 2 , 
Vi = Xj (constant) on dDj, j = 1,2, 



da = ifj^i, 



(3.1) 



I «i(x) - 5 A [l](x) = 0(|x|^ d ) as |x| -> oo. 
Proof. We first observe that 



da = 0, j = 1, 2. 



Since ^1 — IK* is invertible on Hq, there exists a unique solution ) T G %o such 

that 



I-K* 



9 



' S Di [l]\ 



d 



(3.2) 



Let 5ij be the Kronecker symbol and define 



+ 



8i2 



(3.3) 



Since 



we have 



Define 



1, 



-I + K" 



k2 



«i(x) i^Jc^KxHSz^Kx), xGM d \( J D 1 U J D 2 ) 



r 5 



(2)i 



0. 



(3.4) 
(3.5) 



We now show that Vi is the desired function. Because of (|3.4p . we have 







(5 Dl [^ (1) ]+5 D2 [^ j ]) =0 ondD j} j = l,2, 



and hence, + 5d 2 [v 9 !^] is constant in D\ and L>2- Thus, is constant on dD\ 

and <9Z?2- 



Since (-0^ ) T G we have 



(2)i 



d<7 = 0, J = 1,2. 



On the other hand, we have 
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and 

/ da = I 



if j ^ i. Here \dDj\ denotes the area (or length) of dDj. Thus we have the third line in 
(|3.ip . This completes the proof. □ 
As an immediate consequence we have the following theorem. 

Theorem 3.2 The eigenvalue \ o/K* has multiplicity 2. 

Proof. The identity (|3.4p shows that (fj := (tpj ,ipj) T ', j = 1,2, are two eigenfunctions 
of K* corresponding to 1/2. We have from (|3.3p that 

(pfda = \dDj\8ij. (3.6) 

This implies that ip± and ip 2 are linearly independent in TL. Since Hq has codimension 2 
m % and |I - K* is invertible in TLq as mentioned before, the multiplicity of 1/2 is 2. □ 
Using exactly the same arguments one can generalize Theorem 13.21 to the case when 
there are N simply connected inclusions. 

Theorem 3.3 If there are N simply connected mutually disjoint inclusions, then the 
eigenvalue 1/2 o/K* has multiplicity N. 



4 Characterization of the gradient blow-up 



Let <fj, j = 1,2, be the eigenfunctions of K* corresponding to 1/2 introduced in the proof 
of Theorem 13.21 Because of (|3.6|) , if we define g by 



then we have 



Define 



\dDi\ r \dD 2 \ 

[ g^da = l, [ gMda = -l. 
JdDi JdD 2 



?(x) := 5 Dl [ 9 (1) ](x) +S D2 [gM](x), x G R d \ (D 1 UD 2 ) 
Then, q is the solution to 



Ag = in R d \ (Di UD 2 ), 
q = constant on dDi, i = 1, 2, 



Jaoi dJ 1 ) \+ da — 1 ' IdD 2 dJ- 2 ) \+ da ~ ^ 
g(x) = Odxl 1 "^) as |x| ->• oo, 



and 



cM 2 ) 



(4.1) 

(4.2) 
(4.3) 



(4.4) 



(4.5) 
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In fact, since K*[g] = hg, we have 



^(S Dl [g^)+S D2 [ 9 ^})\ 
^_(S Dl [gW}( X )+S D2 [g^})\ 



-V 
2 



7 (2) 



0. 



Thus we have the second line in (|4.4p and ()4.5p . The third line in ()4.4p follows from ()4.5|) . 
Because of (|4.2p . we have 

g(x) = / T(x - y)gW (y) dcr(y) + / T(x - y)g^ (y) da(y) 

= [ [r(x - y) - r(x - yo )] ff « ( y ) My) 

JdD x 

+ / [r(x-y)-r( X -y )y 2) (y) da(y) 
JdD 2 

for any fixed yo- Since 

ircx-tf-rcx-yojl^cixi 1 -" 



as x 



we have the last line in (|4,4D . 

It is known that if D\ and D2 are disks, then the singular function q is given by 
and completely characterizes the blow-up of Vu (see |13j). We have the following theorem 
as a generalization. Here and throughout this paper u\qd denotes the (constant) value of 



u on dDj. 



Proposition 4.1 Assume that D\ and D2 are simply connected domains in M. d , d = 2,3, 
with C 1,a boundaries for some a > 0. The solution u to hi. 5)) can be written as 



u(x) = c e g(x) + 6(x), xel rf \(DiUD 2 ), 
where the constant c t is given by 

u\ 9Dl - u\ 9D2 _ Y,j=i IdDj h ~$3)\+ da 



QldDi - Q\dD 2 



QldDi - Q\dD 2 



(4.6) 



(4.7) 



and 



||V6||L« ( n) < C (4.8) 

for some C independent of e. Here Vt = Vt\\ (D\ U D2) and Oi is an open set containing 
DiUD 2 . 

Remark 4.2 The constant c e in 7| ) may depend on e := dist(Di, D2), but is bounded 
independently of e, and hence the singular function q determines the blow-up ofVu, if D\ 
and D2 are strictly convex and have C 2,a smooth boundaries. This fact will be proved in 
the next section. 



Proof of Proposition 14.11 Let u be the solution to (|1.5p and let 

u := , h : = 



U\dDi 
u\dD 2 



h\dD! 

h\dD 2 
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Since u is constant on dD\ and dD2, it follows from the third line of (|4.4p that 

u\dD! - u\ 9 d 2 = [ u 7rm\+ d<J + I u 7rm\+ d(T = ( u >§)- 

JdD 1 dvW+ J dD2 <9zX 2 J' + 
The representation (|2.4|) implies that 

u = h + S[<p] 

and </9 £ T~Lq. Therefore we have 

(u, g) = (h + S[<p],g) = (h, g) + (<p, S[g]) = (h, g). 



So we have the second identity in f|4.Tj) . 
Let 



6(x) = u(x) — c e g(x) = u(x) r^2-(/(x). 

Then we have 

b\dD! = b\dD 2 - 

So, one can show following the same lines of the proof of Theorem 2.1 in [13] that (|4.8p 
holds. This completes the proof. □ 
It will be shown that the gradient of the singular function q defined as a single layer 
potential of an eigenfunction g of K* blows up as e = dist(Di, D2) — > 0. We now show that 
another eigenfunction does not contribute to the blow-up. Let f = (f^\ f^) T £ T~L be an 
eigenfunction of IK* orthogonal to g with respect to the inner product (|2.12p . namely, 

(S[f],g>=0. 

Then (|4.2|) shows that A^ 1 ) = A^ 2 ^ where is the j'-th component of the (constant) vector 
S[f]. It implies that the function v defined by 

«(x) :=5 Dl [/W](x)+5 D2 [/^](x) 

is constant on dDj, j = 1,2, and satisfies 

v\dDi = v \dD 2 - 
So, \Vv\ stays bounded regardless of e. 



5 Estimates of the potential difference 

We assume for the rest of this paper that D\ and D2 are strictly convex domains in M 2 with 
C 2 ' a boundaries for some a > 0. Let z\ and Z2 be points on dD\ and dD2, respectively, 
such that 

|zi-z 2 |=dist(Di,D 2 ) = e. (5.1) 
We prove the following proposition. 
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Proposition 5.1 Let u be the solution to \ 1.5\) . c t be the constant defined by \4. 7| ), and 

Kj be the curvature of dDj at Zj for j = 1,2. Then c e is bounded independently of e and 



U\dD 1 - U\dD 2 



V«i + «2\/i(l + O a ), 



where 



O, 



0{e a ' 2 ) if a E(0,1), 

0(|Velne|) if a = 1. 



(5.2) 



(5.3) 



We prove Proposition 15.11 after a sequence of lemmas. 

Let Bj be the osculating disk to Dj at Zj so that its radius rj equals to 1/ttj. Let qs 
be the singular function associated with B\ and B2, i.e., the solution to (|4.4|) with Z?i and 
L>2 replaced with i?i and B2. Recall that qs is given explicitly by 



<7b( x ) = ^ (M x 



Pi I 



In |x 



P2| 



(5.4) 



where pi G i?i and P2 G -B2 are the unique fixed points of repeated reflections R\R 2 and 
R2R1, respectively. We emphasize that qs is harmonic in M 2 \ {pi,P2}- 

For the rest of this paper we assume that z\ = (— e/2,0) and Z2 = (e/2,0) after 
translation and rotation if necessary, so that the centers of B\ and B2 are on the real axis. 
In this case, pi and P2 are of the forms pj = (pj,0), j = 1, 2, and it is proved in [23] that 



Pj 



(-r/V2 



T\Ti 
r\ + r 2 



J 



1,2. 



It is also proved using (|5.5p that 



<?b|0Bi - lB\dB 2 



1 



v^vr 



n + r 2 

nr 2 



yfe + 0{e) as e -> 0. 



(5.5) 



(5.6) 



Using (|5.5p one can see that (|1 . 16[) holds. 

Lemma 5.2 There is a constant C independent of e such that 
dq 



01/0') 



(x) 



< C| Vg B (x)| /or all x G 0£>,-, j = 1, 2. 



(5.7) 



Proof. We only prove (|5.7p for j = 2 since the case for j = 1 can be treated in the exactly 
same way. We first assume that D 2 is a disk so that B2 = D2. Let B[ be a disk of radius r[ 
and let 5" be a disk of radius r" (r^ and are independent of e) such that B[ C D\ C -B" 
and z\ G H dB". Let a/ and </' be the solution to (|4.4p with .Cq replaced with and 
-Eq , respectively. Then the monotonic property |21[ Lemma 2.4] yields 



< 



dq 



(x)< 



Q \dB'{ 



1 \dB 2 



dq' 



q'\dB[ ~ Q'\dB 2 I di;W 



(x), x G dD 2 . 



(5.8) 



Because of (|5.6p , there is a constant C independent of e such that 



Q \dB'{ 



Q \dB 2 



<c, 



(5.9) 
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so we have 



Note that 



dq 



</(x) - 



(x) 



<C|Vg'(x)|, x£dD 2 . 



2tt 



fin |x 



Pil 



ln|x - p' 2 |) , 



(5.10) 



(5.11) 



where and p' 2 are the fixed points of the repeated reflections with respect to dB[ and 
dB 2 . Using (|5.5|) one can show that 



|V</(x)| < C|V?fl(x)|, xe5D 2 



(5.12) 



for some constant C independent of e. So we have (|5.7p provided that D 2 is a disk. 

If D 2 is not a disk, we may use a conformal mapping to make it a disk. In fact, if \& 

is a conformal mapping from K 2 U {00} \ U onto M 2 U {00} \ D 2 where U is the unit disk, 
then i& can be extended up to dU as a C 1 function. Therefore, there are constants C\ and 
C 2 such that 

Ci < |V§(x)| < C 2 forallxedc/. (5.13) 

Let 

f = ^ + a 2 , (5.14) 

where a 2 is defined in (|1.4p . Then ^ is a conformal mapping from R 2 \U onto R 2 \ D 2 and 
satisfies (|5.13p . Moreover, by |23l Appendix] and a combination with a linear fractional 
transformation, we can also assume that there are two disks B[ and B'[ of radii independent 
of e such that B[ C ^(Di) C B'{ and <9££ n dB'{ contains the point on <9# _1 (.Di) which 
is the closest to U. Thus we can apply the same argument as above to q o ^ — c q to obtain 
(|5.7p . where c q = linix^oo q o ^(x). This completes the proof. □ 

Lemma 5.3 There exists a positive 5q (independent of e) such that if x £ <9Z)j and |x — 
Zj\ < 5q, then 



and 



Qb(-x) 
8q 



las,- 1 < C-v/i|x - Zj 



x 



X 



+ e 



For any point x G with |x — z,-| > Sq, 



and 



1b (x) -9b|9Bj| < C\/e 
<9g 



(5.15) 
(5.16) 

(5.17) 
(5.18) 



Here, j = 1,2, and the constants C are independent of e and 5q. 



Proof. Assume that j = 2 without loss of generality. There exists So > (independent 
of e) and functions x 2 , xb '■ [— ^o>^o] — > K such that x 2 (0) = e/2, x' 2 (0) = 0, x_b(0) = e/2, 
x' B (0) = 0, and dD 2 and 8B 2 are graphs of £2 and xb for |y| < 5q, i.e., (x 2 (y),y) G <9Z?2 
and {xB(y),y) £ 9-E>2- We then have 



|^2(2/)| < C|2/| 2 and \x 2 (y) - x B (y)\ < C\y 



2+ot 



(5.19) 
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for some constant C since D 2 and B 2 are osculating at z 2 - Since the fixed points pi and 
P2 of the repeated reflections are on the x-axis, we may write pj = (pj,0), j = 1,2, and 
(ELSD holds. 

If \y\ < <5o> then 



i(x 2 (y),y) 
1 



2vr 



ln|(x 2 (y) — pi, 2/)| - ln|(x 2 (y) -P2,y)|) 



1 

2^ 



In |(3;b(7/) -pi,y)| - ln\(x B (y) -p 2l y)\). 



If < \y\ < e 2 ( 2+a > < 5q, there exists x* between x 2 (y) and x B {y) such that 
3{x 2 (y),y) - qs\ 



i&b 2 



<C|x 2 (y)-x B (y)| 



(»* - Pi) 



(x* -p 2 ) 



(x* -pi) 2 + y 2 (x*-p 2 ) 2 + y 2 
< C|x 2 (y) - x B (y)\\pi ~P2\-p2 



ir 



1 



<Cly|^(|y|^« + v^^V% 



If 



where the second to last inequality follows from (|5.5p . 



If e 2(2+ Q ) < |y| < j 0j there exists between p\ and p 2 such that 
QB{x2(y),y) ~ Qb\ 



< C\pi - p 2 \ 



\8B 2 

(x 2 (y) -p*) 



(^s(y) -p* 



(x 2 (y) -p*) 2 + y 2 (x B (y) - p*) 2 + y 2 



< C\pi -p 2 \\x 2 (y) - x B (y)\ 



1 

'y 2 



<C^(|y| 2+a + ^)4<^l2/l a ' 

where the the second to last inequality holds because of (|5.19D . 
If |(x,y) — (e/2, 0)| > 60, one can easily see from ()5.5[) that 



\q B (x,y)\ < C\pi -p 2 \ < Cy/e, 



and hence we have 



i{x 2 {y),y) ~ qB 



8B 2 



< 



Now we estimate 



_8q_ 



'x 2 {y),y)\ + 

on dD 2 . By (|5.7|) . we have 
dq 



qB 



dB 2 



< Cy/e. 



< C\Vq B \ on dD 2 . 



14 



Suppose that \y\ < 5q. Then, we have 



Oq 



< c 

< c 



dx 



(x2(y),y) 



(x 2 (y) - pi) 



(x 2 (y) -p 2 ) 



(x 2 (y) -pi) 2 + y 2 (x 2 (y) - p 2 ) 2 + y 2 
If \y\ < y/e, then \x 2 (y) — Pj\ > Cy/e for j = 1, 2, and thus we obtain 

\Vq(x 2 (y),y)\ <C^. 

If y/e < \y\ < 5o, then it follows that 

\Vq(x 2 (y),y)\<C\ Pl -p 2 \^<C^. 

yZ yZ 

For (x,y) with \ (x,y) — (e/2,0)| > 5q, we have 

\Vq B (x,y)\ <CV~e, 
and (|5.18p follows. This completes the proof. 
Lemma 5.4 We have 



\dD-i - q\dD 2 



1 /j-nr/:,/ 0(e (a+1)/2 ) if a e (0,1), 
— G(|elne|) ifa = h 



as e — > 0. 

Proof. We prove that 



+ 



( e («+i)/2) if a G (0,1), 
0(|elne|) if a = 1, 



as e -> 0. Then (I5T201) follows from (l5T6|) . 
Let 

v(x) := g(x) - Qb(x). 

Since 



dqs 



da 



dqs 



da, i = 1,2, 



the function v satisfies 



Av = mR 2 \(DiUD 2 ), 

u(x) + (?_b(x) — gslaBj = constant on dDi 



dv 



da = 0, i = 1,2, 



v(x) = 0(|x| _1 ) as |x| — > oo. 
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Then, we have 

2 

q\dD 1 - q\dD 2 - (qB\dB 1 - qB\dB 2 ) = - X^ -1 )* ( u ( x ) + qB ^> ~ ^laB^ 



We then obtain from the third line in (|4.4p and the second line in (|5.23|) that 

2 



\dDx ~ q\dD 2 ~ {qB\dB 1 - qB\dB 2 ) =y~] [ (v + q B - qB\dBi) ^~ 



dq 

(0 



do - . 



An integration by parts and the third line in (|5.23p yield 



Let 



q\dDi - q\dD 2 - (tolsBi 



/ (qB-qsldBi) 

JdDi 



dq 



2 /■ 

<=1 79 A 



I \ d( l 



tier. 



+ 



(<7b - ^elasi' 



<9g 



do" 



'|x— Zi|<5 J|x-Zi|>i5o 

Using (|5.17p and (|5.18p we have 

|/ 2 | < Ce. 

To estimate ii, let iV be the smallest integer such that 5q < 2 N yfe. We then have from 
(I5TT5D and ([BTTfil) that 



I'll < 



iV 



< C 



< 



|x- Zi |<v^ - =1 ^2J- 1 v ^<|x-z 4 |<2J v ^ 

/ |x - Zj|dcr + eV / 

■/|x-z,K,/i ^2i"! 



B - </B|a& 



0q 



01/(0 



da 



1 J2J- 1 v ^<|x-z i |<2i v /i 



X — Z< 



|2-a 



^(e + e 1 ^) < 2C a e¥ if a G (0,1), 



C(e + eN) < Celn 



1 



if a = 1. 



□ 



This completes the proof. 

Proof of Proposition 15.11 We first prove that c e is bounded independently of e. For that 
we prove that 

2 

i(h )g )i = y I 

We still assume that zi = (— e/2,0) and Z2 = (e/2,0). Pick a point (c, 0) £ -D2 where 
c is independent of e, and let 



da 



< Cyfe. 



(5.24) 



c 2 y 



(x — c) 2 + y 2 
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(5.25) 



Then, tp is harmonic except at (c, 0) and ip(x,y) = 0(\(x,y)\ 1 ) as |(x,y)| — > oo. Since q 
is constant on dDj, j = 1,2, we have by the divergence theorem that 

^ f* Q f dip ^ f dip 

gL, *^ d " = g L, ^ qia = g* D ' L, ^ da= °- (5 - 26) 

Moreover, one can easily see that there is a constant C > such that 

W(x,y) -y\ < C{x 2 + y 2 ) for all (x, y) G M 2 \ U D 2 ). (5.27) 
Therefore, we have from Taylor's theorem that 



dh 

h(x,y) - h(0,Q) - — (0,0)V(x,y) 



<C(|x|+y 2 



(5.28) 



for all (s,y) G <9L>i U dD 2 . 

Because of the third line in (|4.4p and (|5.26p , we have 

gL,^i + *=gL,( fc -^«')->«w) 



Let 



<9L>, 



0q 



da 



dq 



+ 

|x-zi|<<5 J|x-zi|><5 



da. 



h + h- 



It follows from (f5TT6j) and (pT28j) that 

V<\My)\ + y 2 ) 



\h\ < C 
By (j5TT8|) . we have 

Therefore, we have 



y 2 + e 



da < C 



\ y \<5 V z + e 



^da<C"Ve- 



\h\ < Cyft. 



h - hM - 1(0,0)*)^ 



f/(T 



Similarly, we can show that 



<9£> 2 



„- ft (0,0)-|(0,0)^ 



da 



Hence we obtain (|5.24p . We now infer from (|5.6p and Lemma 15.41 that c e is bounded 
regardless of e. 
Since 



u\dD 1 - u\qd 2 = c € (q\dD 1 - q\dD 2 ) 
by (pr?]) , (|5T2j) follows from Lemma EH 



□ 



17 



6 Estimates of the gradient- Proof of Theorem 11.11 

Proposition 15.11 and (|5.24p show that 

v^(h,g) 



y/e(ni + k 2 ) 



(l + Oo 



where 



Or 



(6.1) 
(6.2) 



0(e a / 2 ) if a 6 (0,1), 

0{\y/e\ne\) if a = 1. 

So, Theorem 11.11 is an immediate consequence of Proposition 14. II and the following propo- 
sition. 



Proposition 6.1 We have 

q(x) = o £ g B (x) + v(x), 



xeR 2 \(DiUD 2 ), 



where 



and 



gjagi — q\dD2 

QBldB! - qB\dB 2 

\\Vv\ 



1 + 



0(e Q / 2 ) 
0(|v^lne|; 

< c 



i/ a 6(0,1) 
ifa=l, 



L^(R 2 \(D 1 UD 2 )) 



(6.3) 

(6.4) 
(6.5) 



/or some constant C independent of e. 



We first fix notation. We suppose that zi = (— e/2,0) and Z2 = (e/2,0) as before. 
There exists do > (independent of e) and functions xi, £2 : [— d~o,d~o] — > K such that 
xi(0) = — e/2, 2^(0) = 0, £2(0) = e/2, ^(O) = 0, and <9Z?i and 3£>2 are graphs of x\ and 
x 2 for |y| < z.e., (xi(y),y) € <9Di and (x 2 (y),y) 6 <9I?2- Since Di and L>2 are strictly 
convex, xi is strictly concave and x 2 is strictly convex. For 5 < 5q, let 

n 5 := {(x, y) e M 2 \ {D 1 U L> 2 ) | x x (y) < x < x 2 (y), \y\ < 6}. 

To prove Proposition 16. 11 we need the following result whose proof will be given in the 
last section. 



Lemma 6.2 If v is a bounded harmonic function in M 2 \ (D x U D 2 ) satisfying 

d 2 v 



Or 2 



[xi{y),y) 



d 2 v 



Or 2 



<M\y\ a ~ l for\y\<5 , 
< M 



(6.6) 
(6.7) 



L°°((dD 1 UdD 2 )\dn So ) 

for some constant M (independent of e), and 

d d 
v(zi) = t^( z i) = v(z 2 ) = g-v(z 2 ) = 0, (6.8) 

where is the tangential derivative on dDi, then there exists a constant C independent 
of e > such that 



\Vv\ 



L°°(R 2 \(DiU£> 2 )) 



< c. 



(6.9) 
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We also need the following lemma. 
Lemma 6.3 There exists a positive constant M independent of e such that 

d 2 q B 



Or 2 



and 



fori = 1,2. 



< M\y\ a ~ l for \y\ < 8 , 



8t 2 



< M 



L°°(dDi\dn 6o ) 



(6.10) 
(6.11) 



Proof. We prove (|6.1U|) and (|6.11|) for i = 2. We use the same notation as in the proof 
of Lemma 15.31 dB<i are given by (x B (y),y) for \y\ < 5q. Let x 2 (y) = (x 2 (y),y) and 
= (x B (y),y). 
Note that 

if \y\ < <5o- Straightforward computations yield 
(<?s(x 2 (y)) - q B (x B {y))) 



1 

2^ 



£(-u 



i+1 



8=1 



(x 2 (y) - • x' 2 (y) (x B (y) - p*) • x' B (y) 



|x 2 (y) - p 



■ 12 



|xs(y) - Pi 



and 



d 2 

^2 (9B(x 2 (y)) - q B (x B (y))) 



W(y)l 5 



i=l 



|x' B (y)|< 



|x 2 (y)-pi| 2 |x jB (y)-p i | 2 



J_ f (x 2 (y) - Pi) -x^y) _ (x B (y) - p,) -x^(y) \ 

+ 2tt^ 1 } V |x 2 (y)- P ,| 2 |xB(y)- Pi | 2 7 



+ T(-ir 

7T ^ — ' 

1=1 

:=/i+/ 2 + / 3 . 



, / ((x 2 (y) - Pi ) ■ x 2 (y)) 2 ((xfld/) - Pi) • y! B {y))< 



|x 2 (y) - Pi| 4 



|x B (y) -Pi| 4 



To estimate 1%, I 2 and ^3, we make some preliminary computations. Since B2 and Z) 2 
are osculating at z 2 , we have 

\x 2 (y)-x B (y)\<C\y\ 2+a . (6.12) 
Since x 2 (0) = x B (0) = e/2 and x' 2 (0) = x' B (0) = 0, we have 
\x2{y)\ + \x B (y)\ <C(y 2 + e), \x' 2 (y)\ + \x' B (y)\<C\y\, \4(y)\ + \x B (y)\ < C, (6.13) 
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and 

\x 2 (y)-Pi\ <C\y\ 2 + V~e, \x B (y)- Pl \<C\y\ 2 + V~t, 3 = 1,2. (6.14) 
It is worth mentioning that the constant C may differ at each appearance. We also have 



\x' 2 {y) ~ x' B {y)\ < Cy 



l+a 



and 



|x 2 (y)- Pl | 2 >C{y 2 + e), \x B (y) - Pi \ 2 > C(y 2 + e), j = 1,2. 
To estimate I±, we write 

\A(y)\ 2 W B (v)\ 2 



X2(y) - Pi| 2 
< \A{y)\ 2 



\xb(v) ~ Pi\ 2 
1 



1 



|x2(y)-Pi| 2 |xs(y)-pi| 2 
Using (JSH2D, dEUD and (fpTTUD we get 



+ 



|x^(y)| 2 -|x^(y)| s 
|x B (y) -Pi| 2 



1 



1 



|x 2 (y)-Pi| 2 |xB(y)-Pi| 2 



We then use (|6.15|) to arrive at 



< 



\x2(y) - x B (y)\\x 2 (y) +x B (y) -2 Pi 
|x 2 (y) - Pi| 2 |x B (y) - pi| 2 



^ c V~e\y\ 2+a + \y\ 4+a <c i 



y 4 + e 2 



\y\ 



l-a ' 



|x 2 (y)| 2 -|x' B (y)| 2 ] < \x'M 2 -x' B {yf 



- Pi 



Thus we have 



|x B (y) - Pi| 2 



< < C 



y 2 + e 



|y| 



l-a ' 



It follows from (|6.17|) that 

(x 2 (y) - Pi) • x' 2 '(y) (x B (?/) - Pi) • x B (y) 



< 



|x 2 (y) - Pi| 2 
x 2(y) • x 2 (y) 



|x 2 (y) - Pi| 2 



|xfl(y) - Pi| 2 



- Pil 2 



+ 

< C 
and hence 



Pi-(xg(y)-x^(y)) 



|x 2 (y) - Pi| 2 



+ 



\*B(y) - Pi| 2 |x 2 (y)-pi| 2 



y 2 + e , y 2 + e , v%| a , ^v%| 2+a + l2/| 4+ ° 



+ 



+ 



2T + e + e y z + e 



+ Ve- 



y 4 + e 2 



<C|y| 



a-l 



\h\ < C\y\ 



a-l 



To estimate ^3, we first write 



((x 2 (y) - pi) • x' 2 (y)) 2 ((x B (y) - Pi) • x' B (y)) s 



|x 2 (y) — Pil 4 



- Pi| 4 



< 



(x 2 (y) - pi) • x' 2 (y) (x B (y) - pi) • x' B (y) 



+ 



|x 2 (y)-Pi| 2 |x B (y)-pi| 2 
(x 2 (y) - Pi) • x' 2 (y) (x B (y) - pi) • x B (y) 



|x 2 (y) - Pi| 2 



|xfl(y) - Pi| 



(6.15) 
(6.16) 



(6.17) 
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One can easily see from (|6.16p that 



x 2 (y) ~Pi)- x' 2 (y) (x B (y) - Pi) • x' B (y) 



+ 



|x 2 (y)~Pi| 2 |x B (y)-pj| 2 
Note that 

(x 2 (y) - pO • x' 2 (y) (xs(y) - p*) ■ x' s (y) 



< 



C\y\ _ 
y 2 + e' 



(6.18) 



|x 2 (y) - Pi 



|xs(y) - Pi | 2 



(x 2 (y) - Pi)x' 2 (y) (^B(y) - PiV B {y) 



|x 2 (y) - Pi | 



|xs(y) 



.•|2 



+ 



x 2 (y)-Pi| 2 |x B (y)-pi 



X2{y)x' 2 {y) x B {y)x' B {y) Pi(x' 2 {y) - x' B (y)) 



|x 2 (y)-Pi| 2 |x jB (y)-Pi| 2 
+ {Pix' B {y) - y) 



|x 2 (y) - Pi| 2 
1 

|xfi(y)-Pi| 2 |x 2 (y)-pi| 2 
We estimate each term using (|6,12p - (|6.16p to have 
(x 2 (y) - Pi) • x' 2 (y) (x s (y) - Pi) • x' B (y) 



|x 2 (y) 



|xfi(y) 



,.12 



< c 



\y\{\y\ 2 + t) , l?/l(|y| 2 + e) , ^|y| 1+Q 



y 2 + e 



y 2 + e 



+ 



y + e 



+ (^ + i)|y| 



^| y | 2+Q + | y | 2+a(y2 + e) 



y 4 + e 2 



<C|y| a . 

Combining this estimates with (|6. 18|) we obtain 

\h\<C\y\ a -\ 

If x € dD 2 satisfies |x — (e/2, 0)| > 6q, it can be easily seen that 

2 v 



8t 2 



(x) 



< M. 



The proof is complete. □ 

Proof of Proposition 16.11 Note that the second identity in (j6.4[) follows from Lemma [57 
and it shows in particular that a e is bounded regardless of e. 
Let 

u(x) = q(x) - a e q B (x), 

and 

u;(x) := — (-y(x) - u(z 2 )). 
Then one can see from the definition fj6.4[) of a e and ()6.19p that 

io(x) = - <7b( x ), x e <9A, i = 1, 2. 

Since Z?i and i?j are osculating at Zj, we have in particular 



(6.19) 
(6.20) 

(6.21) 



It follows from Lemma 16.31 and Lemma 16.21 that 

\\Vw\ 



dm dw 

u>(zi) = -^r( z i) = w(z 2 ) = -q^( z 2) 



< c. 



0. 



lL°° (R 2 \(r>iLLD 2 )) 

Since a e is bounded, we obtain (|6.5p . This completes the proof. 



□ 
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7 The insulating case 



In this section we deal with the case when the inclusions are insulating, namely, the 
problem (|1.20p . We closely follow the argument provided in [T3] . 

Let h 1 - be a harmonic conjugate of h, i.e., h + ih 1 - is analytic. Let u 1 - be the solution 
to (|1.5p with h 1 - in place of h. Then the solution u to (|1.20p is a harmonic conjugate of 
u 1 - in R 2 \ D\ U Z?2 • By Theorem 11.14 we have 

n ± (x) = -^tfi=M=/3 e q B (x) +r(x), x G R 2 \ (D 1 U D 2 ), (7.1) 

yJe{Ki + K 2 ) 

where (3 e is a constant of the form 

Let arg : M 2 \ {(0, 0)} — >■ [— tt, vr) be the argument function with a branch cut along the 
negative real axis, where x = (x\,X2) is identified with x\ + ix 2 . Define 

£b( x ) = ^( ar §( x ~ P 1 ) ~ ar s( x ~ P2) - arg(x- ci) + arg(x - c 2 )), (7.3) 

where Cj is the center of Bj, j = 1, 2. Note that is a harmonic function well defined in 
M 2 \ {B\ U B2) since the jump discontinuity of the argument function across the branch 
cut is cancelled out owing to the fact Pj, Cj G Bj, j = 1,2. Since arg(x — pi) — arg(x — p 2 ) 
is a harmonic conjugate of qs except on the branch cut and |V(arg(x — ci) — arg(x — c 2 ))| 
is bounded independently of e, we arrive at the following result. 

Theorem 7.1 Let u be the solution to il.20\) . Under the same hypothesis as in Theorem 
li.il we have 

u(x) = -^li^^(x)+r(x) > xen\(D 1 UD 2 ), (7.4) 

y/e(Ki + K 2 ) 

where f3 e is a constant of the form ( 7.2 ) and 

l|Vr|| L cc ( n) < C (7.5) 

for some C independent of e. 



8 Boundary value problems 

Let fi be a bounded domain in M? with C 2 boundary. Suppose that O contains two 
perfectly conducting inclusions Dj, j = 1,2, which have C 2,a boundaries for some a > 0. 
We assume that the inclusions are away from dVt, namely, there is a constant cq such that 

dist(Dj,dn) > c , j = 1, 2. (8.1) 
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In this section we consider the following boundary value problem: 



Au = mtt\D 1 UD 2 , 
u = constant on dDj, j = 1, 2, 

du (8 ' 2) 
— — ds = 0, j = 1, 2, 

with either Dirichlet or Neumann boundary conditions on d£l. 

Define a harmonic function h in Q by (|1.22j) . This h plays the role of h in the free 
space problem (|1.5p and we obtain the following result using exactly the same arguments 
as those in Theorem ll.il 

Theorem 8.1 Let u be the solution to lltf.ty) and let h be the function defined by 
Then, under the same hypothesis as in Theorem M.li we have 

u(x) = y^ h ' S \ a e g B (x) + r(x), x £ fi \ (D x U D 2 ), (8.3) 

y'e(/ci + « 2 J 

where a 6 is a constant of the form 
and 

l|Vr|| L =c ( n) < C (8.5) 

/or some C independent of e. 

One can also obtain a similar result for the insulating boundary value problem. 



9 Proof of Lemma 6.2 



We use the same notation as in Section [6j 

Let X2( x ) be a smooth function on dD 2 such that 

X2(x2(y),y) = 1 for \y\ < ^, 

< X2(x 2 (y),y) < 1 for |< \y\ < Sq, (9.1) 
X 2 (x) = otherwise. 

As before we denote the tangential derivative on dDj by J^. Let g 2 + (x) and g 2 _ ( x ) 
be non-negative functions defined for x = (x 2 (y),y) G <9-D 2 > |y| < ooj such that 

92+(z 2 ) = ^±(z 2 ) = <fe_(z 2 ) = ^(z 2 ) = (9.2) 

and 

^ +(x)=mOT /^,o\, ^ (x)=max /_^), l. (9.3) 



<9r 2 <9r 2 <9r 2 <9r 
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Then, g 2+ and g 2 - satisfy 

92+(x 2 {y),y) - g 2 -(x 2 {y),y) = v(x 2 (y),y), \y\ < S . 



(9.4) 



Let V2+, V 2 -, V 2 + and V 2 - be bounded harmonic functions in M 2 \ D 2 which satisfy 
the following Dirichlet boundary conditions on dD 2 : 



V 2+ = X292+ 

V 2 - = X292- 

V 2+ = (1 - X2)max{t>,0} 
V 2 - = - (1 - X2) mm{v, 0} 



on dD 2 . 



(9.5) 



Then by the maximum principle, V 2 +, V 2 -, V 2 + and V 2 - are non-negative and satisfy 



V 2+ - V 2 - + V 2+ - V 2 - = v on dD 2 . 

Let v 2 +, v 2 -, v 2 + and v 2 _ be bounded harmonic functions in 
satisfy the following Dirichlet conditions on dD\ and dD 2 : 



(9.6) 



v 2+ = v 2 - 



v 2+ = v 2 . 



and 




ondDi, 



on dD 2 . 



I 2 \ (Di U D 2 ) which 
(9.7) 

(9.8) 



Since < v 2 ± < V 2 ± and < v 2 ± < V 2 ± on dD\ and <9-D2 ? we have 

< w 2 ± < V 2 ± and < v 2 ± < V 2 ± in M 2 \ (D 1 U L> 2 )- 

Let 

w 2 := ^2+ - «2- + v 2+ - v 2 ^ . 
Then w 2 is a bounded harmonic function in M? \ (D\ U D 2 ) and is such that 

w 2 = on 9Z?i, i«2 = u 011 (9Z?2- 



(9.9) 
(9.10) 

(9.11) 



In the same way, non-negative bounded harmonic functions Vi+, ^i-, ^1+ and V\- in 
i 2 \ {D\ U D 2 ) can be defined so that w± := v\ + — v\- + v± + — v±- satisfies 



w\ = v on dD\, w\ = on dD 2 . 
Then, we have from the uniqueness of the Dirichlet problem that 

v = wi + w 2 mR 2 \ (Di U D 2 ). 
We first estimate | Vv 2+ \ . Thanks to (j6.6[) and (j6.7j) , we have 
d 2 V 2A 



Or 2 



{x 2 (y),y) 



<C\y 



a-l 



for \y\ < 6 , 



(9.12) 



(9.13) 



(9.14) 
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and 



d 2 V 2 . 



dr 2 



< C. (9.15) 

L°°(dD 2 \dn So ) 

Since V 2+ (z 2 ) = £v 2+ (z 2 ) = 0, we have ||V2+||ci,a(aD 2 ) < C- A standard regularity 
estimate for harmonic functions yields, in particular, 



ll^a+llc 



< C. 



(9.16) 



If (x,y) G ILj , then 



\(x,y)-(x(y),y)\ <C{y 2 + e). 
Thus we obtain from (|9.16p and the mean value theorem that 

0< V 2+ (x,y) < C(y 2 + e), (x,y)eU So . 

It then follows from (|9.9|) that 

0<v 2+ (x,y) <C(y 2 + e), (x,y) eIL So . 



(9.17) 



(9.18) 



Let d(x) := dist(x, dD 2 ) and B r (x) be the disk of radius r centered at x. Since v 2 + = 
on dD\ and dD\ is C 2,a , by a standard elliptic regularity estimate we have 



|Vl>2+(x)| < 



C|b2+||L°°(B d(x) (x)n(R 2 \(DiUD 2 ))) 

d(x) 



(9.19) 



If x = (a;, y) G ILj / 2 and x < 0, then there are ci and C2 such that 

ci(y 2 + e) < d(x) < c 2 (y 2 + e). 
If (x', y') G -B rf ( x )(x) n (M 2 \ (Di U D 2 )), then < |y| + c 2 (y 2 + e), thus we have 

\v 2+ (x',y')\<C(y' 2 + e)<C'(y 2 + e). 
We then get from (l9J9]l that 

|vu 2+ (x)| < a 

If x = (x, y) G and x > 0, then we can apply the same argument using the fact 

that V 2 + — v 2+ = on dD 2 to obtain 



|V(V 2+ -is+)(x)| <C. 



We then obtain using (|9.16p that 



|Vl7 2 +(x)| < a 



So, we have 



|v«2+(x)| < c, xen Jo/ , 



(9.20) 



We now estimate |V«2+|- Since || He 1 '" (SZ? 2 ) — we first obtain from the maximum 
principle that 

II^2+||l°° (R 2\ D2 ) < C (9.21) 
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for some C independent of e. Since 

V 2+ (x 2 (y),y) = ifM<V2, (9-22) 

we have 

|W 2+ (x)| < C (9.23) 

for all x G M 2 \ D 2 satisfying 

dist(x,dD 2 \dU So/2 )>S /4. (9.24) 

In particular, we have 

sup |W2+(x)| < C. (9.25) 
xen i()/ 4 

It follows from fljEZg) and ([OS]) that 

0<V 2+ (x,y)<C(y 2 + e), (x,y)GU 5o/4 , 

and from (j9.9|) that 

0<v 2+ {x,y) <C(y 2 + e), (x,y) en jo/4 . 

Since ?52+(xi(y), y) = v 2+ (x 2 (y),y) = if \y\ < So/2, we may apply the same argument as 
for v 2+ to obtain 

|W 2+ (x)|<C, xeII V8 . (9.26) 
In exactly the same way, one can show that 

|V«a_(x)| + |V« 2 _(x)| <C, x G n 5o/8 . (9.27) 

Therefore, we have 

sup |Vu/ 2 (x)| < C. (9.28) 
x en 5o/8 

If x G dD\ \ dllg / 8 , then cZ(x) > C for some C independent of e. Since u> 2 = on (9L>i 
and w 2 is bounded, we obtain 

sup |Vu> 2 (x)| < C. (9.29) 
xeaDi\OT V8 

Let 

V = V 2+ - V 2 - + V 2+ - V 2 - inR 2 \D 2 . (9.30) 
Since V — w 2 = on <9-D 2 , it follows that 

\V(V-w 2 )(x)\<C 

for x G <9-D 2 \ 8ILs /8- Since H^Hc 1 .™^^) * s bounded, || VV r ||£oofiR2\r) 2 ) is bounded, so we 
have 

sup |Vu> 2 (x)| < C. (9.31) 
x.edD 2 \dn So/8 
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Inequalities (EOS]) . (ET2U1) and (EOT!) imply that 



sup 



|Vw; 2 (x)| < C. 



xe0((R 2 \(r>iuz> 2 ))\n, o/8 ) 



We then obtain from the maximum principle that 



sup 



|Vw 2 (x)| < C. 



xe(K 2 \(r>iui?2))\n io/8 



Combining this with (|9.'28p . we readily get 



sup 

xeK 2 \(DiUD 2 ) 



|Vw 2 (x)| < C. 



(9.32) 



One can show in exactly the same way (by switching the roles of D\ and D 2 ) that 
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